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We consider a general class of high-volume, fast-moving production-inventory systems based on both lost- 
sales and backorder inventory models. Such systems require a fundamental understanding of the asymptotic 
behavior of key performance measures under various supply strategies, as well as the pre-planning of these 
strategies. Our analysis relies on a thorongh study of the asymptotic behavior of a random walk with 
power drift, which is of independent interest. In addition to providing key insights, our analysis leads to 
approximations of the corresponding optimization problem that yield simple solutions which are close to 
optimal. We also establish an equivalence between the lost-sales and backorder models when both have the 
same penalty cost that becomes large. 


1. Introduction 

The production-inventory problem we study is that of matching supply and demand over a planning 
horizon of N periods, with the provision that a set of supply decisions must be made before the first 
period. For example, the production-inventory planning horizon might represent a “selling season” 
and the decisions to match supply with demand are made before the season starts. This is not 
merely a quantity decision — how much to stock up so as to fulfill demand over the entire season; 
rather, it is more like a capacity planning decision — how to set the supply or production rates 
as a function of time over the horizon, vis-a-vis the demand forecast. We consider this production- 
inventory matching problem within the context of both lost-sales and backorder inventory models. 

More precisely, consider the following general supply and demand matching problem. Suppose 

Di-\ -h Dn denotes the cumulative demand up to period n = 1,2,..., and assume Var(T)„) = 

(7^ for all n. Further suppose the set of supply decisions before the first period takes the form 
E[T)i] -!-••• + E[Z)„] -b Kan°‘, where a and k are decision variables such that a € [0,1] while k 
can be either positive or negative. This assumed form allows a fair amount of flexibility: one can 
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choose to over-supply (positive k), thus creating safety stock, if the lost-sales penalty is high; or 
to under-supply (negative «;), if the production cost or holding cost is high. We are especially 
interested in asymptotic results. When n is large, demand will approach a normal distribution, 
with a standard deviation equal to a^/n. Shall we match supply with demand by setting a = i? 
What is the corresponding lost-sales quantity? What is the impact on lost sales if a takes on other 
values? What is the right value for n; in particular, should it be positive or negative? 

We note that this general supply and demand matching problem arises in a wide variety of 
production-inventory systems. As such, there is another way to interpret the asymptotics associated 
with large values of N in our model. Instead of viewing this as a prolonged planning horizon, we 
can view the latter as having a fixed length, but with the demand being scaled up. In other words, 

the original total demand for the entire season, D, now becomes Di-\ - \- Dn, with the Z)„ being 

independent copies of D. The capacity problem continues to consist of making the set of supply 
decisions for the entire horizon in advance, with the corresponding deliveries received every period 
to satisfy the demand. Indeed, this alternative view is closer to the class of production-inventory 
applications we have in mind — a high-volume, fast-moving supply and demand context, which 
requires an understanding of the asymptotic behavior of key system performance measures under 
various supply strategies, as well the pre-planning of such strategies. 

As a specific business analytics instance of such capacity problems often encountered within the 
context of workforce management for large services providers, we note that any significant work¬ 
force actions aimed at altering a planned resource supply usually require long lead times, whereas 
customer demand can exhibit volatility and uncertainty over relatively short periods of time. Com¬ 
petition within the marketplace and the dynamics within these companies create constant pressure 
for growth; under this pressure, workforce capacity planning routinely examines resource require¬ 
ments when the demand forecast is scaled up in different proportions. These considerations in part 
motivate the current problem setting. Another application of a similar flavor concerns capacity 
provisioning in cloud computing environments where high-frequency demand having considerable 
uncertainty needs to be fulfilled with a predetermined supply of resources. Any changes in these 
supply decisions are costly and should be avoided. 

Now, let Sn denote the net demand at period n = I,..., A^, namely the difference between demand 
and supply (both in terms of the cumulative quantities) up to period n. In the backorder model, 
where unmet demand in any period will be backlogged and supplied by future supplies, backorder 
and inventory are, respectively, the positive and negative parts of S^, as such, the model is quite 
tractable. On the other hand, the lost-sales model, in which unmet demand in any period will 
be lost, is much harder to solve because the cumulative lost-sales quantity, L„, is the running 
maximum of Sn- Ln = maxo<j<„{5j} (with := 0), which depends on ^i, ^ 2 ,..., Sn- When supply 
is linear, i.e., a = I, then Sn is a random walk and L„ is a well-studied object; in particular, E[L„] 
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is well-known and very accessible. For other values of a, however, E[L„] does not appear to have 
explicit expressions. The only possible exception is a = in which case one can approximate Sn 
by a Brownian motion with square-root drift. We derive such an analysis by relating this Brownian 
motion to the hitting time of an Ornstein-Uhlenbeck process, for which the density function can 
be obtained by inverting Laplace transforms. 

Given this level of difficulty, our approach consists of constructing upper- and lower-bounds for 
E[L 7 v], and analyzing the asymptotic behavior of these bounds for large N. Our hndings for such 
asymptotic behaviors are summarized as follows. When k > 0 (the over-supply case), setting a 
within (i, 1] will result in a bounded E[Ljv], whereas a G [0, will cause E[Ljv] to grow on the order 
of Vn. When k < 0 (the under-supply case), E[L 7 v] will grow on the order of iV“ upon setting a 
within (|, 1]; whereas it will grow on the order of y/N when a G [0, |]. 

In terms of minimizing the total cost over both production cost and lost-sales penalty, our results 
imply that setting a = ^ is the only meaningful choice leading to non-trivial results, independent of 
the values of the penalty and holding costs. Hence, what remains is to hnd the best value of k. We 
first obtain an asymptotically optimal solution through a Brownian approximation of the objective 
function for large N. Then, for each finite N, we examine optimal solutions through the upper and 
lower bounds of E[LAr] as surrogates for the original problem. Both bounds lead to simple optimal 
solutions, and our results include performance guarantees for both problems as approximations to 
the original problem. We also show that these results are readily extended to include holding costs. 
In addition, we establish the equivalence between the lost-sales model and the backorder model 
when both have the same penalty cost that goes to infinity. 

1.1. Related Work 

A brief review of the related literature is in order. To the extent that our set of supply decisions 
has to be made before the season starts, our model resembles the newsvendor model. There are 
important distinctions, however. To start with, the newsvendor decision is a stocking quantity 
decision, whereas we are concerned with the supply dynamics over time — all N periods of the 
planning horizon. Consequently, the lost-sales quantity in the newsvendor model is a single number 
that can be accounted for at the end of the horizon. In our model, lost sales build up over the 
horizon, as demand and supply are realized and evolve. Granted, under our cost model which 
minimizes the total production cost and lost-sales penalty over the horizon, the lost-sales quantity 
E[Ljv] is associated with the end of the horizon, as in the newsvendor model. However, we consider 
inventory holding costs as well the backorder alternative to lost sales, features that are not present 
in the newsvendor model. 

The lost sales problem is known to have a much more complicated structure than the problem 
of backorder, and hence it is more difficult to analyze. Regarding the optimal policy, the exact 
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description is not known even in the simplest case. Karlin and Scarf [12] demonstrate that the 
base-stock policies are not optimal even for systems with a lead time of one period. Monotonicity 
properties of the optimal ordering policy are derived by Morton m, Zipkin [23j and most recently 
Huh [m, each through different methods. Huh et al. uni show that, asymptotically as the ratio of 
unit penalty cost and unit holding cost grows to infinity, the difference in performance between a 
particular base-stock policy and the optimal policy will eventually vanish. Thus, base-stock policies 
can be regarded as asymptotically optimal when there is such a significant imbalance between the 
penalty cost and holding cost. 

Levi et al. [Hj study production and lost sales over multiple periods focusing on the performance 
of a simple heuristic policy that balances the inventory holding cost and lost-sales penalty in each 
period, with the aim of establishing its robustness relative to the optimal policy in a very general 
context. The setting of Levi et al. is dynamic programming: in each period a different decision 
can be made and actions taken accordingly. In contrast, we allow only a single set of decisions 
at the beginning of the horizon; and all subsequent costs due to production, lost-sales or backlog, 
and holding inventory are consequences of this set of decisions. The way we quantify our heuristic 
solution (from the upper-bound problem), however, appears to have a similar flavor to the proven 
“2-approximation” status of the dual-balancing rule in m- Specifically, the performance of our 
heuristic, via the ratio of the upper- and lower-bound solutions, is “2-plus” — where this slightly 
worse performance guarantee than a 2-approximation algorithm can be attributed to the fact that 
there is no recourse over the planning horizon in our model, which, to make matters worse, can be 
infinitely long. A complementary set of results for the lost-sales model has been recently obtained 
by Goldberg et al. |8]. 

Another body of research related to our analysis concerns process flexibility; see the recent paper 
of Chou et al. [7] and the references therein. Specihcally, process flexibility can be modeled by a 
random walk with functional drift, called a generalized random walk in [Tj. When proper functional 
forms are taken, it can be seen that the cost structure under full flexibility has the same form as the 
cost structure of our backorder model, while the cost structure under a so-called chaining option 
corresponds to the structure of our lost-sales model. Upon applying the methodology developed 
in this paper, the effectiveness of chaining can be demonstrated and quantified for more general 
processes. 

1.2. Paper Organization 

The remainder of the paper is organized as follows. In Section [21 we start with a formulation of 
our models — while focusing on the lost-sales model, we also cover the backorder model, followed 
by a derivation of the main performance measures associated with these models. We then present 
preliminaries regarding the normal distribution function and the loss function, and derive estimates 
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for E[Ljv] in the case of a = 1. In Section [31 we continue with estimates of E[Ljv] for all other 
cases, through constructing upper- and lower-bounds and also deriving a corresponding Brownian 
approximation for the case of a = 1/2. Optimization models that minimize total cost over pro¬ 
duction, lost sales (or backorder) and inventory are studied in Section 01 followed by concluding 
remarks in Section [5l 

2. Problem Formulation and Preliminaries 

2.1. Key Performance Measures 

Consider a planning horizon that consists of N periods, indexed by n = 1,..., A^. For each period n, 
let Dn denote the demand, a random variable; and the production or supply quantity, a decision 
variable. At the beginning of the planning horizon, we need to determine the production quantities 
for all N periods, as part of the capacity planning problem. Once this set of decisions is made, 
then in each period n there will be x„ units available along with any units left over from previous 
periods to supply demand £)„. In the lost-sales model, any demand surpassing the available supply 
is lost by the end of each period; whereas in the backorder model, any unfilled demand will be 
backlogged, to be supplied in a later period. Correspondingly, the quantities of interest include the 
following. 

• Ln- the cumulative demand shortfall, in the lost-sales model, up to period n; 

• Bn', the backlogged demand, in the backorder model, at the end of each period n; 

• Hn and H'^: the on-hand inventory left at the end of each period n, respectively for the lost- 
sales and backorder models. 

Proposition 1. Define Sq := 0 and := Di — Xi -I-h — x„. Further define x+ := max{x, 0} 

and x~ := — min{x, 0}. Then, for n = 1,..., A^, we have 

B^ = S-, K = S+, (I) 

Ayi max , Bji (2) 

0<j<n 

Proof. In the backorder model, any unmet demand (if > 0) or any leftover inventory (if Sn < 0) 
at the end of each period n will be carried over; hence, leading to the expressions in ([T|). 

The Hn expression in ([ 2 |) follows from 

n / ^ 

i=l \i=l 

i.e., total production minus total supplied demand equals inventory. 

For the expression in ([ 2 ]), note that the result holds trivially for n = I: Li = max{0, Di — Xi}. 
Suppose it holds for n. Then, we have L„ = for some k <n, which means that there will be no 
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more shortage for periods k + 1 through n. Therefore, the units available to supply D^+i consist 
of Xn+i and the leftover from periods k + 1 through n, namely 

~{Dk+l ~ Xk+l + • • • + Dn — Xn). 

(Note that the above expression is non-negative, as the quantity in parentheses cannot exceed 0 
due to no shortage from periods k + 1 through n.) Hence, at the end of period n-|- 1, shortage will 
occur if 

k^n+l [^n+1 T * * * T Dn ^n)] ^ 0? 

in which case the cumulative shortage will be adding to the above expression, with the sum 
being equal to 5„+i. Otherwise, i.e., if the left-hand side above is < 0, then the cumulative shortage 
stays at S^- In summary, we have 

Ln+i = maxjS'fc, S'n+i} = max{L„, 5„+i} = max {5d, 

0<j<n+l 

which is as desired. □ 

In what follows, we shall focus on since relates directly to L„; whereas and are 
simply the positive and negative parts of the random walk Sn, and thus both are easily accessible 
- refer to m below. 

Since the study of large N asymptotics is our primary objective, we focus on the normal distri¬ 
bution for the demand model. Specifically, let 

/in “t“ 

where are independent and identically distributed (i.i.d.) standard normal random variables, 
and /i„ and a are the mean and standard deviation of D„, respectively. The common a can be 
interpreted as a forecast error; hence, it is independent of the periods. Alternatively, we can allow 
a period-dependent an, as long as the sum ~ he-) it grows at a rate of ^/n. In this 

case, the above expression will still serve as a good approximation due to the central limit theorem. 
On the production side, we write, for n = 1,..., A^, 

Xi-I-|-x„ = /ii-|-h//„-I-Kcrn“, (3) 

where k and a G [0,1] are two policy parameters (or decision variables). Note that the (cumulative) 
production quantity in ([3]) consists of two parts: the first part matches the mean of the (cumulative) 
demand; and the second part can be interpreted as safety stock to offset demand variability. The 
parameter a is often referred to as the “safety factor”. This factor addresses the issue of how the 
safety stock should match up with the demand variability, where the latter grows on the order 
^/n. Recall that k can take on positive or negative values: when it is positive, K,an°' is truly the 
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safety stock; when k is negative, this part is a calculated under-production, suitable for settings 
where the production cost is much higher relative to the lost-sales penalty. Since ([3]) implies that 
the production quantity for each period n is given by 

a:n = /.in + KO-[n“-(n-l)“], n = l,...,A^, (4) 

then the cumulative shortage process L„ can be estimated in different ways depending on the value 
of a. 


2.2. Useful Facts about the Standard Normal Variable 


We summarize below some useful expressions and estimates associated with the standard normal 
variable. Let Z denote the standard normal variate with density and distribution functions denoted 
by (l){x) and 4>(x), respectively. Define the “shortfall function” (or “loss function”) as 


Gix):=E[iZ 



x)4>{z)dz = 4>{x) 


x4)(x), 


(5) 


where $(x) := 1 — <h(x) and the last equality follows directly from 4)'{z) = —z(l>{z). 

Now, we list properties of G{x) that will be useful later for our purposes. 

1. G{x) is decreasing and convex in x (since {Z — x)^ is decreasing and convex in x). 

2. The part of G{x) that has significant curvature is limited to the neighborhood of the origin 
X = 0. For large x > 0, we have G{x) « 0 and G{—x) x. 

3. A direct derivation yields: 

2 f G{x)dx = ^b)-<i>{a) + bGib)-aG{a). (6) 


4. 

0<G(x)<^^, Vx>0. 


(7) 


5. 


x^G(x) 

4>{x) 


To verify dTj) and ([8]), observe that 


as X —>■ -|-oo. 


1 - 


1 \ 1 


x^ 


4>(x) < 4>(x) < —4>(x), 

X X 


X > 0, 


where the upper bound follows from 

X 


( 8 ) 

(9) 


pOO pOO -1 
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and the lower bound follows from 





X 


/2 


The (second) inequality in ([7]) follows immediately from the lower bound in ([2|). Combining the 
upper and lower bounds in dH), we have $(x) ss large (positive) x. This, along with 

rHopital’s rule, helps to verify ([8]). 


2.3. Estimating E[Ln] for a = l 

When a = 1 we have, from dH), = /i„ + na; thus, = —ku + aZ^. In this case, the partial 

sums {^n} constitute a random walk, and E[L„] follows from standard results concerning the 
maximum of a random walk. Specifically, by Spitzer’s Identity (see [20]; also, Ross [IB], Proposition 
7.1.5): 

^ 1 

n—1 

where Sn= — nan + a^/nZ. Hence, 

^[S^]= a^/nE[{Z — Ky/n)~'~] = a^/nG{^^^/n). (11) 

Approximating the summation by integration, we have 

^ -I pN 1 pN 1 Q /'/^\/7V 

J]-E[5+]«/ -E[St]dt = a G{nVi)dt = — G{u)du. (12) 

^ Jo ^ Jo yt ^ Jo 

Therefore, combining the above with dSj), we obtain 

E[LN]^^(^<^>{n^/N)-^'J+aVNG{n^/N). (13) 

When K > 0, the above is increasing in N. As N ^ oo, the second term on the right-hand side 
vanishes, and the hrst term converges to When «: < 0, the second term on the right-hand side 
of ffTB)l quickly becomes a\K\N as N grows, and the first term approaches a constant These 
results then can be summarized in the following proposition. 

Proposition 2. For a = 1 and large N, we have 

E[L 7 v] PS —, when k>0; E[L 7 v] ps —— -+a\K\N, when k < 0. (14) 

ZAv .Z 

Hence, setting a = 1, the expected shortage is bounded by a constant, or grows linearly in N, 
corresponding to a positive or negative n, respectively. 
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3. Estimating E[Ln] for ae [0,1) 

When a < 1, it is generally a difficult problem to calculate the exact value of E[Ljv], and thus 
we develop bounds to reveal its asymptotic order as N ^ oo. We shall initially focus on «; > 0 in 
the next two subsections, and then summarize the corresponding results for k < 0 in the following 
subsection. Lastly, we consider a Brownian approximation for the case when a = |. 

3.1. Case of a G (|, 1) and k > 0 

From ([2|), taking into account Sq = 0, we have 

L„ = max I SA = max 15,^), n = 1,... ,N. 

0<j<n l<j<n ^ 

Jensen’s inequality implies 


E[L^] > max {E[5+]}, 

l<n<N 

where, similar to m, we obtain 

E[5+] = a^/nE.[{Z — = ay/nG{Krf"~^^A ■ 

From d?]), we know that 

a-l/2^ ^ ay/n(j){u) 

ay/nG{Kn°‘ "'A < V ^^ := T , where 

^2^2a-l y2 


(15) 


(16) 


u:=Kn°‘ 


Hence, E[S'+] will decrease to 0 (since (l){u) decreases to 0 exponentially fast as u —)• +oo). The values 
of n that achieve the maximum in the lower bound of ()15p must satisfy the following optimality 
equation (whose solution need not be unique): 


This simplifies to 


G{y) = {2a-l)y^{y) or 0(y) = 2ay$(y), (17) 

and from Q we know that (p{y) ~ y^{y) when y is large. Since a > |, the solution to the above 
equation exists. Indeed, from ([^, we have 


y‘^{y)> ( 1 -^ 

y2 


y>0. 


If we replace y^{y) in (I17p by this lower bound, then the resulting solution 

1 


y:= 


2a-1 


(18) 
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is an upper bound to the solution of the equations in (fT7|) . since </>(y) is decreasing in y for y > 0. 
Therefore, in this case we obtain 

E[LAr] >^max^{E[S'+]}«cj(^0 G(y), (19) 

where y is the solution to (fT7j) and we know y E (0,y). 

For an upper bound, we have = maxo^n^ArlSn} < '^^=1 > which, along with (fTB)) . leads to 

E[Tjv] < a'y^^y/nG{Kn°‘~^) ^ a / y/tG{Kt°‘~^)dt. (20) 

■^0 

Using a transformation of variable u = , together with G{x) < (j){x)/x^ from ([7]), we obtain 


fAT 

'/iG{Kt°‘~^)dt = 


,.Ar““2 


< 


2q; — 1 Jo 
2 

( 2 a - 1)^2 
2 


G{Ku)du 




6-6a 

<p[KUjdU 


3 / 6—6q: . . 

—K ^“-1 / ^20,-1 (pi^vjdv 

1 2o 


2a — 
2 


< --Av cj)(v)dv:= Ga < 00 . 

2a - 1 Jo 

Since > 0, as a E (|, 1), we can summarize our results for this case as follows. 

Proposition 3. For 1/2 < a < 1 and k > 0, we have 
(a) the lower bound: 

N l/(2a-l) 

K 


( 21 ) 


/y* \ l/(2a-l) 

:[Pv]>a(^) G(y:), for iV > y„ := (2a - l)-l/^ 


( 22 ) 


where y* > 0 is the solution to the equation (/)(y) = 2ay$(y) (and hence, independent of N)] 


and 


(b) the upper bound: E[LAr] < aCa, with Ga as specified in (ITT]) . 

Hence, in this case E[Ljv] is both upper- and lower-bounded by constants. 


3.2. Case of a E [0, |] and k > 0 

For a E [0, i], as n ^ 00 , we have G(Kn““^/^) —)• G{0) or G{k) depending on whether a < ^ or 
a = i. Hence, from (fT^ and (fTUD . we obtain, for sufficiently large N, 

E[TAr] > o'G(k)\/]V, whena = -; E[L 7 v] > crG(0)'\/]V, whenaE[0,-). (23) 

For an upper bound, observe that the concave function f{n) = n“ is bounded from below by a 
linear function: 


> —r^n := 0n, for 0 <n< N. 


(24) 
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Defining Wn := —naQn + a know that L^r < max„< 7 v W^- Meanwhile, max„< 7 v Wn can 
be calculated through Spitzer’s identity. More specifically, we have 

E\Wn]= ay/nG{KQ\/n), 0 < re < A^, 


Therefore, making use of ([1|), we obtain 

N 

E[LN]<^-E[W+]^-^[^iKeVN)--'^ +aVNG{KeVN). (25) 

n—1 

As N ^ oo, and according to the definition of 0 in we have Q^/N — 0 if a G [0, \) and 
&y/N —>-1 if a = |. The second term on the right-hand side of (f^ is of order \/]V, whereas the 
first term can be written as — ^]. Since limj,^o(d>(y) — ^)ly = 0(0), we know this 

term is on the order of y/N as well. Hence, 

E[Ln] < a^/N[(f){Ke^/N) + G(re0ViV)] < a^/N[(|){0) + G(0)] = 2cj 0(O)ViV. (26) 

These results allow us to conclude as follows. 

Proposition 4. In the case of a € [0, i] and re > 0, E[Ljv] grows on the order of y/N. In particular, 
we have the lower- and upper-bounds for E[Ljv] in ((251) and (12B)) . 

3.3. Case of re < 0 

First consider the lower hou n ds. For a G (|, 1 ) and re < 0, we have 

] = (T^/nG{^^rG~^/‘^) « cj|re|re“, 

since G{—x) ~ x when x > 0 is large. Hence, for N sufficiently large, 

E[LAr] > max {E[5'+]} cr|re|A^“, foraG(^,l). (27) 

When a G [0, |] and re < 0, the lower bound for re > 0 in (|23l) remains valid; namely, 

E[TAr] > crG(re)-\/]V, whena = -; E[Lj\r]>aG{0)'/N, when a G [0,-). (28) 

Next, we consider the upper bounds. Define Yn:= Sn + recire", for re = 0,1,..., A^. Note that T„ 
is a random walk with zero drift and Tq = 0. Then, 

Ln = max jS'nl = max — Kan°‘} 

0<n<N 0<n<N 
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< max {ynl+ max {—Kan^} 

0<n<N 0<n<N 


0<n<iV 


(29) 


where the inequality follows from the subadditivity of the maximum operator, and the last equality 
is due to K < 0. Applying Spitzer’s identity to the random walk Y^, similar to the derivation of 
E[Ln] in (USD for the case of a = 1, we obtain 


N 


max {y„| 

0<n<N^ \ 


= E 


max 

l<n<Af 


fN 


dt 


{Y+} =^^Gi0)^aGi0) -^ = 2aG{0)VN. 

^ Jo Vt 


Hence, combining this with fl29j) . renders 


E[Ljv] <2cjG(0)ViV + |K|cjiV“, a > 0. 


(30) 


Our results for this case then can be summarized in the following proposition. 

Proposition 5. For k < 0, E[Ljv] grows on the order of cj|k|A^“ when a G (|,1); it grows on 
the order of aG{0)VN when a G [0, i); and it grows on the order of aG{k')\fN when a = |. In 
particular, the lower bound of E[L 7 v] follows ([271) and flTHl) . with its upper bound following (I50ll . 


Finally, when k = 0, we have from (12U1) and fl50)l 

E[LAr] = E' max {Pn}] = 2 (tG(0)\/]V, 


L0<n<Af 

which is also consistent with setting k = 0 in the lower bound (1281) and in the upper bound (|30|) . 


3.4. Brownian approximation for case of a = ^ 

The case of a = I will turn out to produce the right trade-off between cost and service, as we will 
see in the next section, thus yielding the optimal order of production. It is therefore appropriate 
to explore the possibility of a more accurate estimation of L„ in this particular case. 

To this end, we study the continuous counterpart of L„, i.e., a Brownian motion with square-root 
drift. More specifically, E[Lf^] can be approximated as 


E[Ln] 


V^E 


sup < B{s) 
.l<s<t I 



(31) 


where B(t) is a standard Brownian motion. Defining 


:= inf > 1: B{t) = B{1) = — 

I VO' I 


for any a; > 0, we then have 


E 

sup jH(s)--^Vs| 

nOC 

= / 

sup < B(s) - > > X 


.l<s<t t V(J j . 

Jo 

_i<s<t [ yu J 
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P [tx < t] dx. (32) 

Now, define a new process Y{u) := i?(e^“)/e“. One can readily verify, as shown in [3j, that Y{u) 
is an Ornstein-Uhlenbeck process with Y (0) = x. Further define 

Ta_b '■= inf{ri > 0 : y(M) > 6, y(0) = a}. 



We note that, for any t>0, 

P[T,r>e^‘] = (33) 

a fact that has been known since [1]. Our task is thus reduced to the evaluation of Ta^h- When 6 = 0, 
an explicit expression for Ta^b has been derived; see, e.g., [nmaE]. More generally when 6/0 (as is 
of interest here), however, only the Laplace transform of the density of Ta^b is known. Using various 
methods, the moments of this distribution also have been calculated; refer to, e.g., [SaEIlElliBlIlT]. 
In theory, the distribution function of f, can be recovered from all of its moments. In practice, the 
distribution function of Ta^ can be effectively approximated by its finite moments up to a certain 
order. Moreover, as shown in [3], the error bound of such an approximation can be computed 
through numerical methods, such as semi-definite programming. 

We therefore propose to use the first K moments j,... of Ta^b where a = —x and 6 = 
Kj in combination with (1311) . (1321) and fl33p . to obtain a more accurate estimation of in the 
case of a = |. The methods and results in [sainKiT] can be used to calculate the first K moments 
When K <3, we can additionally exploit the closed-form tight bounds in [3] to 
further improve our approximation of the distribution function of Ta^b- (Finding the best possible 
bounds for iL > 4 is NP-hard |3].) As a specific instance of our proposed approach, we consider in 
more detail the case of utilizing the first three moments of Ta^b (he., K = 3). Our starting point is 
the set of closed-form expressions for j,, j,, j, provided in m and parameterized by a = —x 

and 6 = Then, letting X denote a generic random variable that follows the probability 

distribution function of we next obtain closed-form tight upper bounds on both P[A < z] 

and P[X > z] in terms of the expressions for k/^i^-x k/^- ^ore specifically, we 

exploit the following propostion adapted from Theorem 3.3 in [3]. 

Proposition 6. For a constant (i > 0 and a nonegative real random variable X = ^ with 

first three moments Mi = M 2 = M^ , ^ and M 3 = M^ , the functionals P[X > 

(1 -I- 5)Mi\ and P[X < (1 — 5)Mi\ can be approximated as 

P[X>{l + 5)M,]^h{Cl,Dl,,5) and P[X <{l- 5)M,]^ 


respectively, where 




1 

1 Dj^ + (l+S)(Cj^f — S) 

1+5 ■ 15lr + (l+C/)(C/-5)’ 


mm 


s>ci„ 
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(5) — 1 


{ci,+5r 


'M 


{Dl + {Cl + 1){CI + 6)){Dl + {Cl + S)l ’ 


D — 


Ml 

Ml Ms - Ml 
Mf ■ 


(5 < 1 , 


Lastly, we combine Proposition E] and the expressions for ^-x k/vs"’ ^-x ™ m 

together with (fHTI) . ([5^ . ([551) to obtain the more accurate approximation of L„ in the case of a = ^, 
as desired. 


4. Cost Minimization and Asymptotically Optimal Solutions 

In this section we consider instances of two cost minimization models for the lost-sale system of 
interest. For the first model, the total cost consists of both the production cost and lost-sale penalty. 
To be more precise, we minimize cY{,n=iXn +pE[Ljv] where c and p are the per-unit production 

cost and lost-sales penalty, respectively. From ([5]), ignoring the constant term c{pi -|-h Pn), we 

have 


min caKN°‘ +pE[Lj^]. (34) 

For the second model, we incorporate the inventory (holding) cost in each period. Specifically, we 
have 


min caKN°‘+pE[L]\f] + hE[HN], (35) 

(y,K 

where h is the per-unit inventory holding cost. Recall from Proposition [1] that the inventory (at 
the end) of period n is given by — Sn, and therefore E[i7jv] = E[Ln] + KaN°‘. 

Our first step is to find the optimal value of a in these optimization problems. For this purpose, 
let us start with a summary of the results we have derived so far concerning the asymptotics of 
the expected loss sales in the limit as N ^oo. 

(i) When a = 1, E[Ln] is bounded by a constant ^ if k > 0; whereas it grows linearly in N if 

K <0. 

(ii) When a G (i, 1) and k > 0, E[Ljv] is bounded from above and below by constants that are 
independent of N (but dependent on a and k); whereas it grows on the order of if k < 0. 

(hi) When a € [0, E[L 7 v] grows on the order of \/]V, regardless of whether k > 0 or «: < 0. 
From (|54l) . if a G (|, 1] and k > 0, then the second term in the objective will be bounded by a 
constant independent of N; hence, for sufficiently large N, the objective is of order N°‘, i.e., the 
same order as the first term. On the other hand, if a G (|, 1] and k < 0, then the first term will 
decrease while the second term will increase, both on the order of hence, if Op, then 

K = —oo minimizes the objective value, whereas if c < p, then k —)■ 0 minimizes the objective value. 
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Therefore, the only non-trivial solution of the optimization problem is to set a < corresponding 
to the case summarized in (hi) above. Here, the first term of the objective function is of order 
N°‘, which is dominated by the second term of order ^/N (except when a = Consequently, the 
objective value will also be of order y/N, as is further confirmed by replacing E[Ljv] with its lower 
bound via Jensen’s inequality (which applies in all cases): E[Li^] > E[maxi<„< 7 v 5'+] > E[5)()]. This 
leads to the following minimization problem, which is a lower bound of the original problem: 

= cz + pE[{aVNZ — (36) 

To allow for any a G [0,1] and any k, we write 

z = aKN°‘ = ct'/NkN °‘~^ := cj\/iV y, 

and then the optimal solution to minimizing cy + pE[{Z — y)~^] is given by 

„.=*-.(p=$-.(r^). (37) 

With y* being a constant, independent of N, we must have that a = |, and hence k = y*. The 
corresponding objective value then can be expressed as 

ay/Nlcy* +pG{y*)] = a'/N[cy* +p(l){y*)-py*^{y*)] = a\/Np(j){y*). (38) 

This confirms our previous statements; namely, the minimal overall cost cannot be lower than 
order y/N. Moreover, we shall henceforth assume that p> c, unless noted otherwise, because if 
p < c then the lower-bound solution in fl37j) already indicates what will happen: make y as small 
(negative) as possible; whereas the objective value is lower bounded by the expression in (1551) . We 
therefore only need to consider the case of a = Furthermore, upon careful examination of the 
objective of the second optimization model, we can reach the same conclusion. 

In the remainder of this section, we start by considering asymptotically optimal solutions of the 
above optimization models using two different approaches to obtain the desired solutions. We then 
turn to incorporate inventory costs in our optimization models. Finally, we consider the asymptotic 
equivalence of lost-sales and backorder models under certain conditions. 


min cz + pE 


N 


a 




— z 


4.1. Asymptotically Optimal Solutions 

With our focus on a = i, we first obtain an asymptotically optimal solution through a Brownian 
approximation of the objective function for large N. We then obtain optimal solutions through an 
analysis of the asymptotic behavior of our upper and lower bounds of E[Ln]. 
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4.1.1. Brownian Approximation Recall that the optimization problem for the production 
loss trade-off, when a is fixed to be has the form: 


min +pE[Ln]- 

K 

Taking a closer look at the objective function, we conclnde 

^[Ln\ 


(39) 


canM/^+pE[Li,] = (^caK + p^^'^ iV^/T 


Moreover, 


E[L 


N 


Ari/2 


= E 




max 


n<N 


= E 


max 

n<N 


X ”=1 

Ari /2 


where Zi are i.i.d. copies of a random variable following the standard normal distribution. This 
leads to our next result. 

Theorem 7. Let p{k) = sup«;^ Bt — ny/t, where is a standard Brown motion. Then, 

ELi 


lim E 

N—¥00 


max ■ 

n<N 


Proof. We first obtain 


max ■ 

n<N 


E Li 

Ari /2 


iVl /2 


max 


■p{k). 




tG{*.n=l,2....,V} 


tw 
= sup 

tG(0,l] 


iVl /2 


E.':," z. - x-VYm 

NU'i 


From the functional central limit theorem, for the summation of random variables (see e.g., m), 
we can conclude that , , _ 

eS^ 4 - KvW 


iVl/2 


Bt-Vt. 


The desired convergence then follows from the continuous mapping theorem. □ 

The above theorem implies that, for any e > 0 and when N is large enough, we have 

E[Tiv] 


Atl/2 


-p{k) 


< € 


uniformly on compact sets in k. Although there is no explicit formula for p{k), computational 
methods that make nse of the related results in [21I1IIZ112I] together with Monte-Carlo simulation 
can be employed; in particular, these methods can be nsed to obtain accurate estimation of the 
value of p(k) and its derivative in order to solve the Brownian version of our optimization problem: 


min +pap{K). (40) 

K 

Let K* denote an optimum of the problem (|4ni) : of course, k* corresponds to a specihc production 
plan. We next show that this production plan is asymptotically optimal, dehned as follows. 
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Definition 8. A production plan is asymptotically optimal if 

r 

inn sup -< 1. 

N—^oo J /V 

Theorem 9. The optimal solution of (HOll . k*, is asymptotically optimal for the problem ([39]). 

Proof. Let k*{N) denote an optimum of problem (1391) for a given N. Then, for any N, we have 

_ caK*N^^'^+ pE[Ln{k*)] 

Jfj caK*{N)N^l^ +pE[Liy{K*{N))] 

_ caK*M/^+pE[LNiK*)] - +pE[LN{K*iN))]) + + pE[Ln{k*{N))]) 

caK*{N)N^/‘^ + pE[Lpf{K*{N))] 

+pE[Ln{k*)] — (ccrK*(A^)A^A -\-pE[Ln{k*{N))]) + 

“ caK,*{N)N''-/‘^ +pE[Ln{k*{N))] 

We further know that the first term can be bounded by a constant, and thus goes to 0 in the limit 
at N ^ oo. The desired result of asymptotic optimality follows. □ 

It is evident that the asymptotically optimal result is not restricted to the production optimiza¬ 
tion problem, as a similar conclusion can be also reached for the inventory optimization problem. 


4.1.2. Lower- and Upper-Bound Surrogate Problems We now turn to consider an alter¬ 
native approach for obtaining optimal solutions through our upper and lower bounds of E[L 7 v]. 
For finite > 0, we have already analyzed the lower bound surrogate problem in (I38p . We next 
seek to minimize the upper bound with respect to the choice of k. From the upper bound for the 
case a = I and k > 0 in terms of (1331) and (I35]l , upon replacing cay/N and pay/N with c and p, 
respectively, we want to solve the following optimization problem: 


min cK + p ' 

K>0 

The optimality condition is then given by 

P 


1 


-$(k) +G(k) 


(41) 


c+ — (k(I){k) + ^{k) - - ] =p^{k) 


which simplifies to 


ck + pG{k,) = 


_P(^ _ 


:\2 


^{k 


(42) 


Denoting the solution of ((331) by and observing > 0 implies that $(k„) < Further denote 
the objective values of the lower- and upper-bound problems in (1351) and (13T]) by and U“, 
respectively. We write ki :=y* following (|37j) . which yields 


-i/p-c 


Then, from (138|l . (13T]l and (133]l . we have 

_ 2[ck.^+pG{ku)] 

CKi+pG{Ki) 


p 


> 2 , 


(43) 


(44) 
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since ni is the minimizer of cy + pG{y). 

To go in the opposite direction, from (f5Hl) . ([1T]1 and (f^ . we can also write the above ratio as 

Observe that ^[| — ^(a;)] is decreasing in x > 0 (and at x = 0, the limit is 0(0), via rHopital’s rule), 
since its derivative is ^[x0(x) + $(x) — |]. It is then straightforward to verify that x0(x) + $(x) — 
I < 0 for X > 0. We therefore have 

and hence 


20 ( 0 ) 

In fact, we also know that ki< which follows from 


(45) 


_p 

Hi 


(I-*W) 





(since is the minimizer of the right-hand side), together with the fact justified above that 
— $(x)] is decreasing in x. When p = 2c, we have ki = 0; refer to (|151) . Then, we will have = 0 
as well; refer to (|12|). These results render, for this case, = 2V^. 

In summary, we have established the following proposition. 


Proposition 10. Let V* denote the optimal value of the original problem in (fMl) . Let and 
denote the objective values of the lower- and upper-bound problems in (lOHll and (|4T]) . respectively. 
Then, 

_ 2[cKu+pG{ku)] _ 2[I - 4>(At„)] 

CKi+pG{Ki) Ku(j){He) 

where Ki = y* is the lower-bound solution following (|37|) . is the upper-bound solution to the 
optimality equation in (115)) . and Ki< Ku- Consequently, we have 


2 


V 


20 ( 0 ) 

0(Kd 


and 


yl yi yu yu 

- < - < - < 

yu - y* - y* - yi 


Based on the foregoing analysis, we propose to use a = | and k = Ku, the solution to the upper- 
bound problem (for k > 0 ), for our approximate solution to the original optimization problem 
in dSH). When k < 0, we can still use the solution to as a heuristic, simply because the upper- 
or lower-bounds in this case will not lead to meaningful solutions. (For instance, the upper-bound 
problem is min„<o(c — p)k, and thus the solution is: k* = 0 if c < p, and k* = —oo if c > p.) 
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All we need is to remove the constraint k > 0 from the minimization problem in (|4ip , and then 
the optimality equation in (|42p still applies (with the right-hand side increasing in k for k < 0). 
However, it is no longer the case that provides an upper bound. In fact, when k < 0, is a 
lower bound of V*: the constructed linear lower bound to the concave function ^/n, along with a 
negative k, implies subtracting less from the production costs, which in turn must result in lower 
lost-sales penalties. Thus, we have V^/V* < 1, and the performance of as an approximation to 
V*, should be comparable to, if not better than, its performance in the case of k > 0. 


4.2. Incorporating Inventory Cost 

Now, we turn to incorporate inventory costs in our analysis. Referring back to (|35p . let h/N be the 
holding cost per unit of inventory that is held for one period, with h> 0 a constant parameter. The 
reason for including N in the denominator can be explained as follows. The application context we 
have in mind is a fast-moving production-inventory system, where the planning horizon is of a fixed 
length, say one unit, and h is the cost to hold one unit of inventory over this entire horizon. By 
dividing the horizon into N segments (periods), we are effectively using N as a scaling parameter 
on the demand. Hence, as previously noted, instead of demand D over the entire planning horizon, 
we have been considering for each period n = with Dn=D and independent of n. 

Another way to motivate this approach is due to the fact that, if h is not scaled by N, then the 
inventory holding cost, because of its cumulative nature over time, will dominate all other costs 
when N grows. 

Therefore, using a = |, the inventory cost is given by 

h ^ 2 

— ^ a^/n[G{K) -I- k] PS -ha'/N[G{K) + k] 


for the lower bound; and similarly, for the upper bound, we have 


h 

N 


N 

'^ay/n 

n—1 



<&(«;) +G{k) + k 


-hay/N (— - — -|-G(k)- t-; 


V K L2 


refer to ([251) and (f25l) . Upon replacing |/i by h, the inventory costs derived above are easily 
incorporated into both the Brownian approximation and the lower- and upper-bound surrogate 
problems of the previous subsections: simply replace c by c -|- h, and p hy p + h. 


4.3. Asymptotic Equivalence of Lost-Sales and Backorder Models 

Next, we consider another type of asymptotic results. Suppose the penalty p is large, say p^ oo, 
whereas N is fixed. We first ask the question of whether it is possible to keep the expected penalty 
cost pE[L]\f] finite (for a given N). 
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To this end, let us consider the lower-bound solution with a = | and k, = Ki = y*] for simplicity, 
we shall write y instead of y*. Recall from (l37)l that y = when inventory cost is included 


into the model, and thus p - 


■ oo is equivalent to y 
P^[Sn] = crVnpG{y) = 


oo. Then, for any n, we have 
cG{y) 


Hy) ’ 


for which applying I’Hopital’s rule yields 


lim SM = lini = lini i ^ 0. 

y-yoo y-^oo y-yoo y 

Namely, ^£[5+] —)• 0 for every re¬ 
in fact, with S^{y) := ay/n{Z — y)+, this strategy also leads to E[5+(y)] = ay/nG{y) —)• 0 as 
y —)■ oo for every re. Because S^{y) is decreasing in y, monotone convergence implies that 5+(y) —)• 0 
as y —)• oo, for every re = 1,..., Hence, we must have Lj^{y) = maxi<„<jv as y-^oo, 

and thus E[Ljv(y)] ^ 0 as y —>■ oo. Since the optimal strategy (to the original problem with inventory 
cost included) can do no worse, we must have 


mm a 


T^fN{c + h)K + (p-E/ i)E[Ljv(k)]^ 

< cr\/]V(c-E/i)y-E (p + /i)E[Ljv(y)] = a\/N{c +h)y + o{y), 


where y = as specified above. 

On the other hand, based on the results established earlier, we know that the right-hand side 
above is also a lower bound of the objective function, because (with k* denoting the minimizer of 
the first expression) 


min (cj\/]V {c + h)n + {p + /i)E[L7v(re)] j 

> a\/N{c + h)K*+ {p + h)E.[Sjj{K*)] 

> ay/N{c + h)y+{p + h)E[Sji{y)] = cr\/]V(c-E/i)y-E o(y), 


where the second inequality follows from y being the minimizer of the lower-bound problem. We 
therefore have 

min ^cr\/iV(c-E/i)K-E (p + /i)E[Ljv(k)]^ = cr-\/]V(c-E/i)y-Eo(y). (46) 

Next, let us consider the backorder model. Recall that = 5"+ is the number of backlogged units 
of demand and that = S~ is the inventory, both at the end of period re; refer to Proposition 
[TJ Similar to the lost-sales model in the previous subsection, let ^ and ^ be the per-unit backlog 
penalty and inventory holding costs in each period, respectively. Then, the (original) optimization 
problem can be expressed as 

/ , w ,, N \ 

mm ck-E —J]]E[5+(k)]-E —J]]E[5'-(k)] J =aVN min[{c +h')K + {b + h')G{K)], (47) 

\ n=l n=l / 
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where |6 and |/i' are replaced by b and h' (analogons to the approach taken with h in the previous 
subsection), respectively, and taking into account S'” = 5+ — = 5+ + ay/uK. Let y' denote the 

optimal solution to ()T7ll . clearly rendering y' = Suppose we set b = p and h' = h, and let 

oo. Then, y' = y and the left-hand side of (|47p has the following asymptotics for large y: 

aVN min (^cK + bE[S:^{K)] + h'[S~{K)]^ = {c +h)y' + o{y') = (c + h)y + o{y). (48) 

Upon comparing (flB)) and we have established the following result. 

Proposition 11. As p^oo, the objective values of the lost-sales and backorder models become 
equivalent in that their ratio approaches 1, provided the lost-sales and backlog penalty costs are 
both equal to p, and the purchasing and holding costs are also respectively equal in the two models. 

5. Conclusions 

A general class of high-volume, fast-moving production-inventory problems are studied. Although 
we focus on the case in which the production is set to be the average demand supplemented 
by a certain amount of safety stock, our assumptions on this safety stock is quite general. In 
the asymptotic regime, we identify the right order of the safety stock by studying random walks 
with power drifts. We also derive bounds and approximations for functionals of key performance 
measures so that the production planning can be optimized under different settings. Finally, our 
analysis provides another means and settings for establishing the asymptotic equivalence of lost 
sale and backorder models, as observed by several authors using different means and settings. 
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